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ABSTRACT 

We present an XMM-Newton observation of the massive edge-on Sb galaxy NGC 2613. We 
discover that this galaxy contains a deeply embedded active nucleus with a 0.3-10 keV lumi- 
nosity of 3.3x10''° ergs s~^ and a line-of-sight absorption column of 1.2 x 10^'^ cm^^. Within 
the 25 mag arcsec"^ optical B-band isophote of the galaxy, we detect an additional 4 sources 
with an accumulated luminosity of 4.3 x 10'^^ ergs s~^. The bulk of the unresolved X-ray emission 
spatially follows the near-infrared (NIR) K-band surface brightness distribution;the luminosity 
ratio Lx/Lk ^ 8 x 10~^ is consistent with that inferred from galactic discrete sources. This 
X-ray- NIR association and the compatibility of the X-ray spectral fit with the expected spectrum 
of a population of discrete sources suggest that low-mass X-ray binaries (LMXBs) are the most 
likely emitters of the unresolved emission in the disk region. The remaining unresolved emission 
is primarily due to extraplanar hot gas. The luminosity of this gas is at least a factor of 10 
less than that predicted by recent simulations of intergalactic gas accretion by such a massive 
galaxy with a circular rotation speed Vc ^ 304 km s"-' (Toft et al. 2002). Instead, we find that 
the extraplanar hot gas most likely represents discrete extensions away from the disk, including 
two "bubble-like" features on either side of the nucleus. These extensions appear to correlate 
with radio continuum emission and, energetically, can be easily explained by outflows from the 
galactic disk. 
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X-ray observations of extraplanar hot gas (T 
> 10® K) around nearby edge-on disk galaxies 
are essential in the study of the galactic ecosys- 
tem in many aspects, particularly the disk-halo 
interaction. Such observations have helped es- 
tablish the prevalence of galactic superwinds in 
starburst galaxies, e.g., NGC 253 (Strickland et 
al. 2000, 2002) and NGC 4666 (Dahlem, Weaver 
& Heckman 1998), among others. Extraplanar 
X-ray-emitting gas has also been detected un- 
ambiguously around several "normal" late-type 
galaxies with little evidence for nuclear star- 
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bursts: NGC 891 (Sb; Bregman & Houck 1997), 
NGC 4631 (Scd; Wang et al. 2001), NGC 3556 (Sc; 
Wang, Chaves & Irwin 2003) and NGC 4634 (Scd; 
Tiillmann et al. 2006). In these galaxies (except 
for NGC 4634 which currently lacks direct evi- 
dence), extraplanar hot gas is clearly linked to 
outflows from recent massive star-forming regions 
in galactic disks. The global X-ray properties of 
extraplanar gas in these "normal" star-forming 
galaxies, when scaled with the star formation rate 
of the host galaxies, appear similar to those found 
in starburst galaxies (Strickland et al. 2004a, b; 
Wang 2005). Nevertheless, this needs to be con- 
firmed by extended X-ray observations of "nor- 
mal" star-forming galaxies. 

On the other hand, current galaxy formation 
models also predict the existence of hot gaseous 
halos surrounding present-day disk galaxies, which 
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arise from gravitational infall from the intergalac- 
tic medium (IGM; e.g., Toft et al. 2002 and refer- 
ences therein). The predicted extraplanar X-ray 
luminosity strongly depends on the mass of the 
host galaxy. X-ray observations thus have long 
been expected to detect such gaseous halos around 
nearby massive, typically earlicr-typc disk galax- 
ies. However, there is so far little direct obser- 
vational evidence for the presence of this kind of 
X-ray-emitting halo. Benson ct al. (2000) ana- 
lyzed X-ray emission from the outer halos (> 5') 
of primarily two early-type spirals NGC 2841 (Sb; 
~ 15 Mpc) and NGC 4594 (Sa; - 25 Mpc), using 
ROSAT PSPC observations. No significant dif- 
fuse emission was detected, although the upper 
limits to the diffuse X-ray luminosities arc consis- 
tent with the current predictions (Toft et al. 2002). 
Therefore, more dedicated searches for the X-ray 
signals of IGM accretion around disk galaxies are 
needed. 

Here we present a study of an XMM-Newton ob- 
servation toward NGC 2613, an edge-on Sb galaxy 
with "normal" star formation. We focus on prob- 
ing the spatial and spectral properties of its large- 
scale X-ray emission. This galaxy (Table 1) is a 
good candidate to probe the presence of hot gas, 
in the sense that: 1) it is very massive and thus 
expcctcid to contain a large amount of hot gas; 2) 
its high inclination (~79°) allows the possibility 
of detecting extraplanar emission, either from a 
halo of accreted gas or a large-scale outflow; 3) 
its moderately large distance (25.9 Mpc) places 
the galaxy and its ~50 kpc vicinity in the field 
of view (FOV) of a typical XMM-Newton observa- 
tion, offering a good opportunity of studying the 
large-scale distribution of gas, and 4) it is known 
to show extraplanar features at other wavebands, 
specifically the radio continuum and HI (Chaves &: 
Irwin 2001; Irwin & Chaves 2003) as well as earlier 
HI observations (Bottema 1989). 

2. Observations and Data Reduction 

We obtained two XMM-Newton observations 
on NGC 2613. The first observation (Obs. ID 
0149160101), taken on April 23/24, 2003 with 
a 40 ks exposure, suffered heavily from back- 
ground flares. Consequently, another observation 
(Obs. ID 0149160201) with a 33 ks exposure was 
taken on May 20/21, 2003. In our analysis, we 



only used data obtained from the second observa- 
tion. 

We used SAS, version 6.1.0, together with the 
latest calibration files for data reduction. For the 
MOS data, we selected only events with patterns 
through 12 and applied flag filters XMMEA_EM, 
XMMEA_2, XMMEA_3 and XMMEA_11. For the 
PN data, we selected only events with patterns 
through 4 and applied flag filters XMMEA_EP, 
XMMEA_2. XMMEA.3 and XMMEA.ll. Ac- 
cording to light curves of the MOS and the PN, 
we further excluded time intervals with high 
background rates by setting good time interval 
thresolds of 2.0 cts/s for the MOS in the 0.3-12 
keV range and 6.0 cts/s for the PN in the 0.3-15 
keV range, respectively. The resulting net expo- 
sure time is 23.5 ks for the MOS and 18.7 ks for 
the PN. We used the skycast program to generate 
the "blank-sky" background for our observation. 
Same event filters were applied to the "blank-sky" 
event files, resulting in exposures of 791, 759 and 
294 ks for the MOSl, M0S2 and PN, respectively 
We then constructed images and exposure maps 
at the 0.2-0.5, 0.5-1, 1-2, 2-4.5 and 4.5-7.5 keV 
bands for each detector. For spectral analysis, we 
selected only events with FLAG = 0. 

3. Analysis and Results 

3.1. Discrete X-ray Sources 

We perform source detection on the PN images 

of the soft (S, 0.5-2 keV), hard (H, 2-7.5 keV) 
and broad (B=S-|-H) bands. As detailed in Wang 
(2004), the source detection procedure, optimized 
to detect point-like sources, uses a combination 
of detection algorithms: wavelet, sliding-box and 
maximum likelihood centroid fitting. The source 
detection uses a detection aperture of the 50% 
PSF encircled energy radius (EER). Multiple de- 
tections with overlapping 2a centroid error circles 
are considered to be the same source, and the cen- 
troid position with the smallest error is adopted. 
The accepted sources all have a local false detec- 
tion probability P < 10"''. 

We detect a total of 67 discrete sources on the 
PN images, 5 of which are located within the 
Ib = 25 mag arcscc"^ isophotc of the galaxy. We 
have also examined the source detection in the 0.2- 
0.5 keV band, but found that the results are con- 
taminated by many artifacts, especially near the 
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CCD boundaries. Using the same detection pro- 
cedure, we detected 46 sources in the combined 
images of MOSl and M0S2, but no new source is 
found within the FOV of the PN. Tabic 2 hsts the 
sources detected from the PN data, the locations 
of which are shown in Fig. 1. 
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Fig. 1. — EPIC-PN intensity image in the 0.5- 
7.5 keV band after flat-fielding. An adaptively 
smoothed background has been subtracted from 
the image to highlight the discrete sources, which 
are marked by circles with radii of twice the 
50% EER. The source numbers (Table 2) are also 
marked. The ellipse illustrates the optical Ib = 
25 mag arcsec'^ isophote (7^2x1^8) of NGC 2613. 
The dashed circle outlines the region where the lo- 
cal background spectrum is extracted for spectral 
analysis. 

The prominent nucleus of NGC 2613 is readily 
seen in Fig. 2a. The nucleus is the only source 
detected within the central 30", suggesting that 
this emission may represent an AGN. To inves- 
tigate this further, we perform a spectral analy- 
sis on the nuclear emission. Due to the relatively 
low spatial resolution and short exposure time of 
the observation as well as the relatively large dis- 
tance to NGC 2613, the spectrum extraction of 
the nucleus is a compromise between having bet- 
ter counting statistics and suffering less contami- 
nation from non-nuclear emission around the nu- 
cleus. To assess the contribution from non-nuclear 
emission, we extract two spectra from circles with 



radii of 10" and 16" around the galactic center, 
for each detector. The 10" (16") radius corre- 
sponds to a physical scale of ~1.2 (2.0) kpc at 
the distance of NGC 2613 and represents an en- 
closed energy fraction (EEF) of ~0.50 (0.65) in 
the PN and ~0.55 (0.70) in the MOS. We extract 
a background spectrum for each detector from a 
circle with a 2. '5 radius at ~ 6' to the south of the 
galactic disk (Fig. 1). This background region is 
chosen because it shows a low X-ray intensity and 
is at an off-axis angular distance comparable to 
that of the disk. We then bin the source spectra 
to achieve a background-subtracted signal-to-noise 
ratio greater than 2. 

k-i' 
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Fig. 2. — EPIC-PN intensity images around the 
galactic center in the 2-7.5 (a) and 0.5-2 (b) keV 
bands. The nucleus is prominent only in the hard 
band. The plus sign marks the optical center of 
NGC 2613. The dashed circle with a 16" radius 
illustrates the region for the spectrum extraction. 

The spectra from the 16" circle show a promi- 
nent broad bump at energies ^ 4-5 keV (Fig. 3), 
which can also be seen in the spectra from the 
10" circle. This is further evidence that the nu- 
clear region contains an AGN since the spectra 
can naturally be explained by a combination of a 
heavily absorbed power-law from an AGN and a 
softer contribution from non-nuclear emission. We 
jointly fit the PN, MOSl and M0S2 spectra in 
the 0.3-12 keV range with XSPEC. Owing to the 
lack of obvious features below 3 keV in the spec- 
tra, a composite model of wabs[zwabs(PL)+PLj 
is applied, where the first power-law component 
(PLl) with intrinsic absorption characterizes the 
nuclear emission and the second power-law com- 
ponent (PL2) represents the non-nuclear emis- 
sion. In the fit, we require that the amount of 
foreground absorption be at least that supplied 
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by the Galactic foreground (as specified in wabs: 
ISIhi > 6.8x1020 cni-2). Fit results for the 10" 
and 16" spectra are consistent with each other 
within the uncertainty ranges. We list in Ta- 
ble 3 the fit results to the 16" spectra, the im- 
plications of which will be discussed in § 4.1. All 
quoted errors in the tables are at the 90% confi- 
dence level. The best-fit two-component model is 
shown in Fig. 3. 

We note that in the PN spectrum there is some 
hint of excess over the best-fit two-component 
model at ~0.9 keV, which might be physically due 
to the presence of diffuse hot gas around the nu- 
cleus. We thus try to include a thermal plasma 
component (APEC in XSPEC) in the fit to probe 
the existence of an additional thermal component. 
While the best-fit temperature of this component 
is, as expected, ~0.9 keV, the fit is not signifi- 
cantly improved according to an F-test, and the 
range of the fitted temperature could not be well 
constrained. Therefore we consider this potential 
thermal component insignificant in the spectra. 
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Fig. 3.— EPIC spectra (Black: PN spectrum; 
Red: MOSl spectrum; Green: M0S2 spectrum) 
of the nucleus of NGC 2613 extracted from a 16" 
circle. The best-fit model (solid curve) consisting 
of two absorbed power-law components (dashed 
and dash-dot curves, respectively) is also plotted. 
The lower panel shows the data-to-model ratios. 

Given the relatively high source detection limit 
(~ 2 X 10^^ ergs s~^) and the relatively low spa- 
tial resolution, the bulk of expected galactic X-ray 
sources of NGC 2613 are still embedded in un- 
resolved emission. Indeed, only four sources (in 
addition to the nucleus) are detected within the 
D25 ellipse of the galaxy. We jointly fit the ac- 



cumulated spectra of these four sources extracted 
from individual detectors. Circles with radii of 
twice the 50% EER around individual sources are 
used to accumulate the spectra. Corresponding 
background spectra are extracted from the source- 
removed D25 ellipse. We find that an absorbed 
power-low model fits the spectra well, giving a 
best- fit photon index of LSOtpJ? and an intrinsic 
0.3-10 keV luminosity of 4.3 x 10^^ ergs 8"^ The 
slope of the power-law is typical for X-ray binaries. 
Given their high luminosities (> lO'^^ ergs s~^), 
on average, these sources are likely ultraluminous 
X-ray binaries. 

3.2. Unresolved X-ray emission 

3.2.1. Spatial distribution 

As an overview. Fig. 4 shows the large scale 
0.5-2 keV X-ray emission around NGC 2613 and 
its similarity with the optical disk. Along the ma- 
jor axis of the disk, the emission appears rather 
smooth and is confined within ~2.'5 from the 
galactic center. Whereas along the direction per- 
pendicular to the disk, some extended features are 
present, forming well-defined structures. To the 
north of the disk (Fig. 4) is a "bubble-like" fea- 
ture, which is referred to as the 'north bubble' in 
the following. This feature follows the minor axis 
fairly well and has a maximum extent of ~ 100" 
(13 kpc) from the nucleus. Immediately on the 
opposite side of the nucleus is another extension 
that is somewhat smaller, reaching ^1'. This fea- 
ture will be called the 'south extension'. A third 
feature protrudes from the south side of the ma- 
jor axis but west of the minor axis. This feature 
will be called the 'south-west feature'. Finally two 
very large extensions are seen to the east of the 
major axis with emission peaks located at RA= 
08''33™37^ DEC= -22°59'19" (north arc), and 
RA= 08''33™33", DEC= -23°0'51" (south arc). 
These have the appearance of arising from the 
eastern tip of the X-ray disk and will be called 
the 'eastern extensions' consisting of northern and 
southern arcs. The above features are labelled in 
Fig. 4 for ease of reference. Most of these fea- 
tures extend significantly beyond the optical disk. 
Given the extent of the X-ray emission along the 
major axis and the inclination of the disk, the pro- 
jected in-disk emission along the minor axis should 
be within 45", whereas most of the extended fea- 
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turcs show an extent larger than 1'. Therefore, we 




Fig. 4.— EPIC-PN 0.5-2 keV intensity contours 
overlaid on the digitized sky-survey (first gener- 
ation) image of NGC 2613. The X-ray intensity 
is adaptively smoothed with CIAO csmooth to 
achieve a signal-to-noise ratio of ~3. The contour 
levels are at 4.1, 5.0, 7.0, 8.8, 13, 16, 22 and 32 
xlO~^ cts s~^ arcmin"^. The dotted ellipse rep- 
resents the optical Ib = 25 mag arcsec"^ isophote 
of the galaxy. The dashed ellipse (5' x 3') illus- 
trates the region where spectra of unresolved X- 
ray emission are extracted. 

We remove the detected discrete sources, ex- 
cept for the nucleus, to study the spatial proper- 
ties of the unresolved soft X-ray emission. The nu- 
cleus is kept because it is heavily obscured (§ 3.1) 
and contributes few photons to the emission be- 
low 2 keV. For each faint source with a count rate 
{CR) < 0.01 cts s~^, we exclude a circular region 
with a radius of twice the 50% enclosed energy ra- 
dius (EER). For sources with CR > 0.01 cts s~^, 
the radius is multiplied by an additional factor 
of l-Mog(Ci?/0.01). The choice of the adopted 
radii is a compromise between removing the bulk 
of the source contribution and preserving a suf- 
ficient field for the study of unresolved emission. 
Thus 75-80% of photons from each detected source 
is removed according to our criteria. 



We construct the "blank-sky" background- 
subtracted, exposure-corrected 0.5-2 keV inten- 
sity profile from the source-removed PN image, as 
a function of off-center distance along the minor 
axis. The full width along the direction parallel 
to the major axis used for averaging the inten- 
sity is 5' 0.7 -D25), approximately the maximal 
extent of the unresolved emission along the ma- 
jor axis (Fig. 4). The profile is shown in Fig. 5. 
We characterize the profile by an exponential law: 
I{R) = lo e~l^l/^°, where \z\ is the vertical dis- 
tance from the center, zq is the scale height and 
Iq is the central intensity. A constant intensity 
is included in the fit to account for the discrep- 
ancy between the local background and the sub- 
tracted "blank-sky" background, and it turns out 
to be negligible. The results are listed in Table 4. 
The best-fit model is also plotted in Fig. 5 (black 
dashed curve) . There is an excess over the best- fit 
model at a vertical distance of ~l'-lf5, which can 
be attributed to the extraplanar features apparent 
in Fig. 4. The total count rate produced by this 
extraplanar excess is ~ 5.2 x 10^"^ cts s^^. 

The unresolved emission should consist of 
two components: 1) emission from truly diffuse 
gas, and 2) collective discrete contributions from 
sources below our detection limit plus some resid- 
ual counts spilled outside our source removal re- 
gions. To constrain the source component, we 
assume that it follows the distribution of the NIR 
K-band light of the galaxy, which can be deter- 
mined from the 2MASS K-band map (Jarrett et 
al. 2003). We remove from the map bright fore- 
ground stars and convolve it with the PSF of 
the PN. Circular regions used for removing the 
discrete X-ray sources are also excluded from the 
map. The K-band vertic;al profile is then produced 
in the same manner as for the X-ray profile. To 
compare the two profiles, we normalize the K-band 
profile to match the X-ray intensity at the center 
(Fig. 5). This requires a normalization factor of 
3.0x10"^ cts s"^ arcmin-2/(MJy sr"^). Within 
a vertical distance of 1', the K-band profile traces 
the X-ray profile quite well. An excess over the 
K-band profile, however, remains at a vertical dis- 
tance of ^l'-lf5. This is evidence for the presence 
of extraplanar X-ray-emitting gas. 

We also construct a 0.5-2 keV radial intensity 
profile for the unresolved emission. Elliptical pho- 
tometry is applied, with a minor-to-major axis 
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Fig. 5.— EPIC-PN 0.5-2 keV intensity profile 
(black crosses) along the direction perpendicular 
to the disk of NGC 2613. A "blank-sky" back- 
ground has been subtracted and discrete sources 
have been removed, except for the nucleus. The 
full width used for averaging the intensity is 5'. 
Spatial binning is adaptively adjusted to achieve 
a count-to-noise ratio greater than 12, with a 
minimum step size of 12". The black dashed 
curve is a fit to the profile using an exponential 
law plus a constant local background. The red 
solid curve shows the 2MASS K-band profile con- 
volved with the PN PSF and normalized by a 
factor of 3.0x10^3 cts s^^ arcmin^Vl^Jy sr^^) 
(see text). A constant of 10~^ cts s~^ arcmin"^ 
has been added to all data points to avoid nega- 
tive values improper for a logarithmic plot. 

ratio of 0.29 and a position angle of 113° (Jar- 
rett et al. 2003). We again fit the radial pro- 
file with an exponential law and list the results 
in Table 4. The X-ray profile together with the 
best-fit model is plotted in Fig. 6. Also plotted 
is the normalized 2MASS K-band radial profile 
produced in the same way as is done for the X- 
ray profile. The K-band and X-ray profiles closely 
trace with each other within a semi-major radius 
of ^ 3'. At semi- major radii ^A'-6', corresonding 
to semi-minor radii '^I'.l-l'.l, a bump is present. 



again most likely due to the extraplanar features. 
The total count rate of this bump is consistent 
with the extraplanar excess seen in the vertical 
intensity profile. 
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Fig. 6. — Similar to Fig. 5 but for the radial pro- 
file, generated with the elliptical photometry (see 
text). A constant of 10^"^ cts s^^ arcmin^^ has 
been added to all data points to avoid negative 
values improper for a logarithmic plot. 

3.2.2. Spectral properties 

Guided by the X-ray morphology (Fig. 4), a 
galactocentric 5' x 3' ellipse with a position an- 
gle of 113° is adopted to extract the spectra of 
unresolved emission for individual detectors. The 
limited number of counts prevents us from further 
dividing the spectra according to different regions 
of interest, e.g., the disk and the halo. Discrete 
sources are removed in the way as described above, 
except for the nucleus, for which we exclude a cir- 
cular region with a radius of 40" (~87% EER) to 
further reduce the contamination from the hard 
nuclear emission. We use the same background 
spectra as applied in § 3.1. 

We jointly fit the PN, MOSl and M0S2 spec- 
tra in the 0.3-8 keV range (Fig. 7). OveraU, the 
spectra are much softer than those from the nu- 
clear region, showing clear features at ^ 0.9 keV, 
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corresponding to the Fe L-shell complex, and at 
~ 0.5 keV. These features further indicate the 
presence of diffuse hot gas. At energies above 2 
keV the spectra are dominated by a collective con- 
tribution from unresolved discrete sources, most 
likely LMXBs (see § 4.2). We account for this 
contribution with a power-law (PL) with a pho- 
ton index fixed at 1.56, as found by Irwin, Athey 
& Bregman (2003) for the accumulated spectra of 
galactic LMXBs. This PL, combined with a ther- 
mal plasma component (APEC) is then used to 
fit the spectra. Foreground absorption is again re- 
quired to be at least the Galactic value. The fit is 
initialized by fixing the plasma metal abundance 
at solar. The model, however, yields a poor fit 
to the spectra, in particular failing to simultane- 
ously account for the features at ~0.9 keV and 
'^O.S keV. By allowing the metal abundance to be 
fitted, we obtain a statistically better fit but the 
resulting abundance is low (< 0.05 solar). Such an 
extremely sub-solar abundance is unphysical and 
practically often encountered in the X-ray spec- 
tral analysis for galaxies (e.g., NGC 253, Strick- 
land et al. 2002; NGC 4631, Wang et al. 2001). 
We thus add to the model a second APEC com- 
ponent. A two-temperature plasma is effective 
in characterizing the diffuse spectra of some star- 
forming disk galaxies (e.g., Strickland et al. 2004a; 
T-iillmann et al. 2006). The abundances for both 
thermal components are fixed at solar. The fit 
is acceptable, resulting in a cool component with 
kT ~0.08 keV and a hot component with kT -0.8 
keV. The 0.3-10 keV X-ray luminosities are 5.2, 2.6 
and 1.1x10'^^ ergs s~^ for the three components, 
namely the discrete sources, the low temperature 
gas and the high temperature gas, respectively. 
The above results are summarized in Table 5. We 
adopt the two-temperature fit as the best-fit model 
in the following and plot it in Fig. 7. 

4. Discussion 

4.1. The nature of the nuclear X-ray emis- 
sion 

Our two-component spectral fit for the nucleus 
(§ 3.1; Table 3) indicates that the intrinsic neu- 
tral hydrogen column density is '^1.2x10^'^ cm~^. 
This is much higher than the beam-averaged HI 
column density of —2 x 10^^ cm^^ found by Chaves 
& Irwin (2001), but it is typical of values found for 
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Fig. 7. — EPIC spectra of unresolved X-ray emis- 
sion (Black: PN spectrum; Red: MOSl spectrum; 
Green: M0S2 spectrum) of NGC 2613 and the 
best-fit wabs(PL+2APEC) model. The spectra 
are adaptively binned to achieve a background- 
subtracted signal-to-noise ratio better than 2. The 
lower panel shows the data-to-model ratios. 

molecular circumnuclear disks. Ott et al. (2001), 
for example, find a molecular column density of or- 
der 10^3 cm^2 foj. 4945. Molecular data are 
not yet available for NGC 2613, but our results 
suggest that a substantial molecular component 
should be present in this galaxy. 

For the nuclear component, the modelled in- 
trinsic flux given in Table 3 leads to an intrin- 
sic X-ray luminosity of —3.3x10'*° ergs s^* in the 
0.3-10 keV range. The photon index of this com- 
ponent is —2, a typical value found in the X- 
ray spectra of AGNs (e.g. Pellegrini, Fabbiano & 
Kim 2003). No radio core was detected by Irwin, 
Saikia & English (2000), putting a 3cr upper limit 
of 4.5 X 10^^ ergs s^* Hz^^ on the radio spectral 
power at 1.425 GHz within the same 16" region. 
Using the above luminosity over the 0.3-10 keV 
range, we derive an upper limit of a = 0.62 on 
the energy spectral index (Si, oc h'~°') between 
the radio and X-ray bands. Although the X-ray 
nucleus is heavily obscured, these values neverthe- 
less suggest that the energy spectral index is likely 
flat or possibly rising at the low frequencies, a fact 
again consistent with the interpretation of the nu- 
clear source as an AGN. Thus, we conclude that 
the nuclear X-ray source represents an AGN in 
this galaxy, the first evidence that this is the case. 

The non-nuclear component, characterized by 
the second power-law (PL2), shows a photon in- 
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dex of ~ 1.7 and an intrinsic 0.3-10 keV lumi- 
nosity of 3.4x1039 ergs s"^. Irwin et al. (2003) 
showed that the accumulated spectra of LMXBs 
in carly-typc galaxies can be uniformly described 
by a power-law model with a best-fit photon in- 
dex of 1.56 ± 0.02. By using a sample of nearby 
galaxies of various morphological types, Gilfanov 
(2004) studied the relation between the collective 
luminosity of LMXBs and the K-band luminos- 
ity, Lk, of the underlying stellar content. He 
found that, for LMXBs with luminosity higher 
than 10^'' ergs s~^, their collective luminosity Lx 
= (3.3- 7.5) X 10^9 ergs s-iLK/lO"L0,if. As- 
suming that the spatial distribution of LMXBs 
follows that of the K-band star light (Jarrett 
et al. 2003), we estimate that the collective 
luminosity of LMXBs within the 16" circle is 
~(1.8 - 4.2)xl03^ ergs s~-^. Thus the non-nuclear 
component is consistent with the collective emis- 
sion of unresolved LMXBs. We note that high- 
mass X-ray binaries (HMXBs) are expected to 
be present in star-forming disk galaxies and their 
composite spectral properties are somewhat sim- 
ilar to that of the LMXBs, thus the collective 
contribution of HMXBs may also be partly re- 
sponsible for the non-nuclear component. We 
show below that in NGC 2613 the relative con- 
tribution of HMXBs is small as compared to that 
of LMXBs. 

4.2. The collective X-ray emission of dis- 
crete sources 

It is known that X-ray binaries, including 
LMXBs and HMXBs, dominate the X-ray source 
populations with luminosities > 10^''' ergs s"-'^ in 
galaxies. Owing to their distinct evolution time- 
scales, the numbers and thus the collective con- 
tributions of long-lived LMXBs and short-lived 
HMXBs to the X-ray emission of a galaxy are ex- 
pected to be proportional to its stellar mass and 
star formation rate (SFR), respectively. Colbert et 
al. (2004) analyzed Chandra observations of X-ray 
sources in a sample of nearby galaxies of various 
morphological types and SFRs. They found that 
the collective X-ray luminosity of point sources 
Lxp is linearly correlated with the total stellar 
mass Mi, and the SFR of the host galaxy as 

Lxp (ergs s"^) = (1.3 ±0.2) x 10^^ M* (Mq) 

-F(0.7 ± 0.2) X 10^^ SFR (Mq yr'^). (1) 



We use this relation to assess the relative im- 
portance of LMXBs and HMXBs in contributing 
to the X-ray emission of NGC 2613. The total 
stellar mass can be estimated from the K-band lu- 
minosity Lk and the B — V color index via (Bell 
& de Jong 2001) 

log(Af*/LA-) = -0.692 + 0.652(B - V), (2) 

where Lk is in units of the K-band Solar lumi- 
nosity. The SFR can be estimated from the far- 
infrared (FIR) luminosity LpiR via (Kennicutt 
1998) 

SFR = 4.5 X 10"^* Lfir (ergs s'^). (3) 

LpiR is measured according to (Lonsdale, Helou 
& Good 1989) 

Lfir = 3.1 x lO^^ (2.58 Sqo + Swo), (4) 

where D is the distance of the galaxy in units of 
Mpc, 5*60 and Sioo are the flux densities in units of 
Jy at 60 /xm and 100 /xm, respectively. With the 
available photometric data for NGC 2613 (Table 
1), we estimate that the total stellar mass is 2.1 x 
10^1 Mq and the SFR is 4.2 Mq yr"!. Based on 
Eq. (1), the contributions of LMXBs and HMXBs 
to the collective X-ray emission of discrete sources 
is ~ 2.7 X 10^° ergs s'^ and ~ 2.9 x 10^^ ergs s-\ 
respectively, with the latter being about 10% of 
the former. 

In the disk of NGC 2613, we find that the 
0.5-2 keV unresolved emission is spatially cor- 
related with the K-band star light. Therefore, 
the normalization factor for the K-band profile 
(§ 3.2.1) should represent the collective X-ray 
emissivity of the underlying old stellar popula- 
tion. Using the power-law model given by Ir- 
win et al. (2003) for the accumulated spectrum 
of LMXBs, we convert the observed 0.5-2 keV 
count rate into the intrinsic luminosity in the 0.3- 
10 keV band). The K-band flux density is also con- 
verted into intrinsic luminosity according to the 
2MASS K-band photometry. The normalization 
factor, 3.0x10-" cts s'^ arcmin-2/(MJy sr^^), is 
then equivalent to an X-ray emissivity of Lx = 
4.2x10^^ ergs s~''^Lk/{10^^Lq^k), or a luminos- 
ity ratio of Lx/Lk ^ 7.5 x IQ—*. Gilfanova 
(2004) found that the collective X-ray luminos- 
ity of galactic LMXBs is related to the under- 
lying K-band luminosity following Lx = (3.3 — 
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7.5) xlO-'''' ergs s-1Lx/(1O"L0,k), i.e., a lumi- 
nosity ratio of Lx/Lk ~ (5.8 - 13.2) x 10"^. 
Therefore the collective X-ray emissivity of unre- 
solved discrete sources inferred for NGC 2613 is 
consistent with that of the galactic LMXBs, and 
we conclude that the collective X-ray emission of 
LMXBs dominates the soft emission of NGC 2613 
in its disk region. 

4.3. The origin of extraplanar gas 

Presence of diffuse gas in NGC 2613 is evident 
by the soft X-ray excess over the K-band light. We 
consider two possible origins of the diffuse gas: 1) 
the continuously accreted IGM (Toft ct al. 2002) 
and 2) the outflow from the galactic disk (Irwin & 
Chaves 2003). 

4-3.1. An accreted gaseous halo? 

Toft et al. (2002) calculated global X-ray prop- 
erties (e.g., luminosity, effective temperature and 

intensity distribution) of hot gaseous halos, based 
on their simulated galaxies. The predicted lu- 
minosity strongly depends on the circular speed 
of the host galaxy. The most massive galaxies 
in their simulations have circular speeds similar 
to that of NGC 2613 (~ 300 km s'^). The pre- 
dicted 0.2-2 keV luminosity for such a galaxy is 
~ 8 X ergs s^^ (Fig. 3 in Toft et al. 2002). 
Prom our best-fit spectral models of the spec- 
tra of unresolved emission, we derive an intrinsic 
0.2-2 keV luminosity of ^8x10^*^ ergs s^-^ for the 
sum of the thermal and power-law components, 
and ~6x 10^^ ergs s^^ for the thermal components 
only. We note that the unresolved emission out- 
side our spectral extraction region contributes lit- 
tle to the total luminosity. The simulated luminos- 
ity of gas emission by Toft et al. (2002) is at least 
an order of magnitude higher than the observed 
value for NGC 2613. We therefore conclude that 
the simulations as presented by Toft et al. (2002) 
substantially over-predict the X-ray emission from 
the cooling inflow of the IGM, if this is what is oc- 
curring in NGC 2613. 

This over-prediction is related to the so-called 

over-cooling problem in current theories of galaxy 
formation. We speculate that the over-cooling 
problem is a result of an inappropriate treatment 
of stellar and/or AGN feedback. For example, the 
mechanical energy input from Type la supernovae 



(SNe) is typically not included in galaxy forma- 
tion simulations, partly due to the difficulty in 
treating the astrophysics related to gaseous flows. 
Qualitatively, Type la SNe, which tend to occur 
in low-density hot environments, provide an es- 
pecially effective mechanism for large-scale dis- 
tributed heating, required to reduce the cooling 
of gas in galactic bulges and halos (Tang & Wang 
2005). Massive stars in galactic disks may also 
serve as sources of mechanical energy that could 
produce outflows into halos and help slow down 
the cooling of the accrected gas. For example, 
with a star formation rate of ^ 4.2 yr~^ 
for stars between 0.1 and 100 Mq, and assum- 
ing a Salpeter IMF and that stars with mass 
> 8 Mq become core-collapse SNe, the rate of to- 
tal energy release from the star-forming regions of 
NGC 2613 is Lsnii ~ 1-0 x 10^^ gj-gs s'^ Our 
spatial and spectral analyses suggest that the ex- 
traplanar gas is responsible for the thermal emis- 
sion (§ 3.2) and has a total 0.3-10 keV luminosity 
of ^ 5 X 10'^^ ergs s~^. Thus, SNe can provide 
enough energy to explain the X-ray emission of 
the extraplanar gas in NGC 2613. 

4.3.2. Multiwavelength extraplanar feat,ures 

In Fig. 8a and b, we compare the X-ray emis- 
sion with the radio continuum emission and HI to- 
tal intensity emission, respectively. Of the two ra- 
dio images, the radio continuum morphology more 
closely resembles the X-ray morphology in the 
sense that: a) the north bubble has a radio contin- 
uum counterpart; b) the south extension also has 
a radio continuum counterpart; c) the south-west 
feature shows a small radio continuum protrusion; 
and d) the peaks of the large eastern extensions 
(north and south) also show radio emission. 

The HI total intensity map shown here does 
not show all of the extended features identified by 
Chaves & Irwin (2001), but two of their features, 
Fl and F2 clearly extend above and below the 
galactic plane and are labelled in Fig. 8b. These 
two features might be related with the northern 
and southern arc of the eastern extensions seen in 
the X-ray. 

It is not wise to read too much into these corre- 
lations, given the limited S/N of the maps. How- 
ever, the relationship with the radio continuum is 
sufficiently strong that the X-ray emission in the 
extraplanar features, representing hot diffuse gas. 
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Fig. 8. — (a) VLA C+D configuration continuum contours overlaid on the same X-ray intensity image (grey 
scale) as contoured in Fig. 4 and in (b) of this figure. The contour levels are 0.18, 0.27, 0.56, 0.84, 1.1, 
1.7, 2.3, 3.2 mJy bearn"^ and the beam is 22"xl5" at a position angle of —8? 2. A few X-ray extraplanar 
features are labelled (see text), (b) The same X-ray intensity contours as in Fig. 4 overlaid on a greyscale 
image of the total intensity VLA C+D configuration HI map. The grey scale range (shown with a square 
root transfer function) is in units of 10^ Jy beam"^ m s~^ and the beam is 47"x32" at a position angle of 
— 8?2. Fl and F2 refer to two HI extensions identified by Chaves & Irwin (2001). 



is very likely associated with the radio continuum 
emission which represents predominantly the non- 
thermal component. 

Wc further consider the energetics of a spe- 
cific feature, namely the 'north bubble', which is 
the only extraplanar feature that can be cleanly 
isolated from the ambient emission. Guided by 
Fig. 9, we approximate the volume occupation of 
the bubble by a cylinder with 1' in diameter and 
0.'8 in height, the center of which is 1'.2 above the 
galactic center. Hence the volume of the bub- 
ble is ~ 2.7 X 10^ kpc^. We find a total 0.5- 
2 keV count rate of 2.4x10""^ cts s~^ within the 
bubble. In the best-fit model to the spectra of 
unresolved emission, the high and low tempera- 
ture components predict a 0.5-2 keV count rate 
of 7.5x10"^ cts s~^ and 2.0x10"^ cts s~^, respec- 
tively. Therefore the north bubble is unlikely to 
be due to the low temperature component alone. 
Instead, it could be dominated by the high tem- 
perature component. Taking an effective temper- 
ature of ^ 0.8 keV, we estimate the mean den- 
sity of the bubble to be ~ ry"^/^ x 10~^ cm~^, 
where r] is the filling factor of the hot gas inside 
the bubble. The total thermal energy of the bub- 
ble is Eth — 3.677"^/^ X 10^^ ergs s~^, and the work 



done to steadily lift up the bubble against gravity 
is Eg ~ 1.7?7~^/^ X 10^^ ergs s^^, given the grav- 
itational potential introduced by the exponential 
disk of the galaxy (Irwin & Chaves 2003). Given 
the morphology and the position of the bubble, 
we speculate that it was produced near the nu- 
clear region, either by a nuclear starburst or the 
AGN. If the bubble's total amount of thermal and 
gravitational energy is obtained from a starburst, 
it takes a time of t ~ {Eih + Eg)/{fLsNii) — 
1.7?7~^/^ X 10^ yr to form the present structure, 
where / is a geometrical factor taken to be 0.1 to 
reflect a fractional star formation rate of the cen- 
tral 1 kpc in the disk. This timescale is typical 
for massive stars to become SN explosions. On 
the other hand, assuming the flux density of the 
AGN follows Si, oc v~°' between the radio and 
X-ray bands, the total bolometric luminosity over 
this frequency range is 4.5 x 10^° ergs s~^ with 
a = 0.62 adopted (§4.1). It is uncertain what frac- 
tion of the AGN energy can be taken to energize 
the ambient gas, but we consider the AGN might 
also be capable of producing this feature. For ex- 
ample, the locations of the north bubble; and south 
extension immediately on either side of the nucleus 
are reminiscent of extraplanar loops or lobes seen 
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in nuclear outflow galaxies like NGC 3079 (e.g. 
Cecil et al. 2002) which is known to have an AGN. 




RIGHT ASCENSION (J2000) 

Fig. 9. Blow-up of the X-ray-emitting bubble to 
the north of the nucleus. The same X-ray intensity 
image as in Fig. 4 is used. Contour levels are at 5, 
5.2, 5.4, 5.6, 5.9, and 6.2 xlQ-^ cts s'^ arcmin'^. 

5. Summary 

We have analyzed an XMM-Newton observa- 
tion of the massive edge-on Sb galaxy NGC 2613. 
We find a deeply embedded AGN in this galaxy. 
The X-ray spectrum of this AGN can be char- 
acterized by a power-law model with a photon- 
index of ~2 and a 0.3-10 keV intrinsic luminosity 
of 3.3x10''° ergs s"!. Linking the X -ray spectral 
properties of the AGN with the current upper limit 
at radio frequencies indicates a spectral flattening 
of the AGN at low frequencies. 

The 0.5-2 keV unresolved X-ray emission is 
found to closely trace the near-IR emission in the 
disk region, and the X-ray to ncar-IR luminosity 
ratio is consistent with that inferred from galactic 
LMXBs. These two facts together indicate that 
the bulk of the unresolved emission is produced 
by the old stellar population of the galaxy, pre- 
dominantly LMXBs. 

A few extraplanar diffuse X-ray features are 
present in addition to the collective emission from 
discrete sources traced by the near-IR light. These 
features can be explained by the presence of hot 



gas, which can be spectrally characterized by a 
two-temperature plasma with kT ~0.08 keV and 
~0.8 keV. The total X-ray luminosity of hot gas is 
at least an order of magnitude lower than that 
predicted by current simulations of IGM accre- 
tion based on disk galaxy formation models. Thus 
the extraplanar features are very unlikely to result 
from IGM accretion. 

Instead, morphologically most of these extra- 
planar features have extended radio counterparts, 
which are believed to arise from disk-related 
events. Also, energetically the extraplanar fea- 
tures can be generated by either supernova explo- 
sions or the AGN, the latter possibly related to 
the bubbles above and below the nucleus. There- 
fore, we conclude that the extraplanar features are 
most likely formed from outflows from the galactic 
disk. 

Our observation suggests that a proper inclu- 
sion of galactic feedback is essential, not only to 
understanding galaxy formation, but also to its 
continued evolution. NGC 2613 and galaxies like 
it provide nearby laboratories that may help to 
understand the over-cooling problem existing in 
current galaxy formation simulations. 

For J. I., this work has been supported by the 
Natural Sciences and Engineering Research Coun- 
cil of Canada. 
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Table 1 

Basic information of NGC 2613 



Parameter NGC 2613 



Morphology" SA(s)b 

Optical size" 7f2 x IfS 

Inclination angle*" 79° 

Position anglc"^ 113° 

Center position" R.A. 08^33™22!84 

(J2000) Dec. -22°58'25'/2 

B-band magnitude " 11.16 

V-band magnitude " 10.25 

K-band magnitude ^ 6.82 

60 /im flux (Jy)'' 7.48 

100 /xm flux (Jy)"^ 25.86 

Circular speed (km s"-^)*" 304 

Distance (Mpc)^ 25.9 

(1' = 7.53kpc) 

Redshift" 0.00559 

Galactic foreground A^hi (10^° cm.-^)" 6.8 



References. — a. NED; b. Chaves & Irwin (2001); c. Jarrett 
et al. (2003); d. Sanders et al. (2003); e. Dickey & Lockman 
(1990). 
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Table 2 
XMM-Newton Source List 
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1.9 


7.04 ± 1.17 




-0.07 ±0.20 


B, S 


42 


J083334.38-230542.1 


2.1 


9.12 ± 1.40 




-0.43 ±0.16 


B, S 


43 


J083336.73-225922.2 


2.9 


2.76 ±0.85 






S 


44 


J083337.40-225124.4 


3.8 


3.27 ±0.99 






s 


45 


J083338.72-230155.1 


2.5 


3.82 ±0.88 






B, S 


46 


J083340.39-230529.7 


3.5 


4.32 ± 1.07 






B, S 


47 


J083340.73-225056.0 


4.0 


4.72 ± 1.18 






B 


48 


J083342.87-225030.0 


5.0 


3.73 ± 1.14 






S 


49 


J083350.48-231150.8 


5.2 


9.28 ±2.19 






B, S 


50 


J083351.30-225656.7 


4.3 


3.31 ± 1.04 






S 
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Table 2 — Continued 



Source 


XMMU Name 


(") 


CR ( cts ks-i) 


HR 


HRl 


Flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


51 


J083352.03-225306.0 


3.0 


10.61 ± 1.63 




-0.33 ±0.15 


B, S 


52 


J083352.44-230516.5 


3.3 


6.49 ± 1.38 






B, S 


53 


J083352.84^231135.6 


5.4 


7.58 ± 2.05 






B, S 


54 


J083353.06-225728.9 


2.5 


7.78 ± 1.32 




-0.16 ±0.18 


B, S 


55 


J083354.10-230211.9 


3.1 


5.35 ± 1.21 






B, S 


56 


J083356.35-230249.5 


3.4 


5.16 ±1.30 






S, B 


57 


J083358.74-230730.8 


4.3 


7.31 ± 1.71 






B, S 


58 


J083359.63-230915.8 


4.0 


25.05 ± 3.05 




-0.06 ±0.12 


B, S 


59 


J083406.25-224808.9 


5.7 


46.66 ± 4.55 


-0.05 ±0.14 


0.38 ±0.11 


B, S, H 


60 


J083410.07-225808.2 


4.4 


6.46 ± 1.55 






B, S 


61 


J083416.92-225508.4 


5.6 


6.17 ± 1.85 






S 


62 


J083417.03-230448.5 


5.8 


7.14 ±2.04 






B 


63 


J083422.96-225632.9 


4.8 


67.11 ±5.02 


-0.11 ±0.11 


0.10 ±0.09 


B, S, H 


64 


J083426.43-225921.1 


5.8 


11.11 ±2.42 






S, B 


65 


J083428.17-225454.8 


6.0 


26.10 ±3.35 




-0.07 ±0.14 


B, S 


66 


J083428.60-225508.9 


7.0 


18.48 ±2.93 




-0.07 ±0.19 


H 


67 


J083430.98-225843.4 


6.3 


11.55 ±2.42 






B, S 



Note.— The definition of the bands: 0.5-1 (SI), 1-2 (S2), 2-4.5 (HI), and 4.5-7.5 kcV (H2). In addition, 
S=S1+S2, H=H1+H2, and B=S+H. Column (1); Gcnerie source number. (2): XMM-Ncwton X-ray 
Observatory (unregistered) source name, following the XMM-Newton naming convention and the lAU 
Recommendation for Nomenclature (http://cdsweb.u-strasbg.fr/iau-spec.html). (3): Position uncertainty 
(Ic) calculated from the maximum likelihood centroiding. (4): On-axis source broad-band count rate — 
the sum of the exposure-corrected count rates in the four bands. (5-6): The hardness ratios defined as 
HR = (H - S2)/(H + S2), and HRl = (S2 - Sl)/S, listed only for values with uncertainties less than 0.2. 
(7): The label "B", "S", or "H" mark the band in which a source is detected with the most accurate 
position that is adopted in Column (3). 



15 



Table 3 

Spectral fit to the nuclear emission" 



Parameter Value 

xVd-o.f. 64.6/65 

''Nhi 6.8 (<11.4) 

"Na/ 12.3t^^7^ 

Photon index (PLl) 2.ltJ j 

Photon index (PL2) l-Vto'g 

%.3-io (PLl) 41 

%.3-w (PL2) 4.2 



Note. — a. See text for model de- 
scription; b. Foreground column den- 
sity in units of 10^" cm~^, mini- 
mum sets at the Galactic foreground 
value of 6.8; c. Intrinsic column den- 
sity in imits of 10^^ cm~^; d. In- 
trinsic 0.3-10 keV flux in units of 
10~^* ergs cm~^ s~^. 



Table 4 

Fit to 0.5-2 keV surface intensity distributions 



Parameter Vertical distribution Radial distribution 



xVd-o.f 


42.4/24 


65.2/23 


"lo (10~^ cts s~^ arcmin"^) 


15.3tli 


57.0t«1 


^zq (ro) (arcmin) 


n qo+0.04 
"••^^-0.04 




''lb (10~^ cts s~^ arcmin"^) 


0.0i};° 


'^O.O 



Note. — a. Central intensity; b. Scale height (length) of the exponential 

law; c. Local background intensity above the already subtracted "blank-sky" 
background; d. Fixed at the value obtained from the vertical fit. 
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Table 5 

Spectral fit to the unresolved emission 



Modd 




<nj " 


a'' 


Touipcraturo'' 


A1)Uii(l;uico 


Flux'' 


\-/d.o.[. 


PL+APEC 

PL+APEC 
PL+2APEC 


6i 
6.8 
6.8 


I (<7.9) 

(<17.9) 
(<2L9) 


1.56 

1.56 
1.56 




1.0^ 
0.02 (<0.05) 
1.0^= 


7.5 (PL), 1.3 
3.2 (PL), 4.4 
6.5 (PL), 3.2, 1.4 


81.1/66 

57.9/65 
62.5/64 



Note. a. Column density in units of 10^" cm"^, minimum sets at the Galactic foreground value of 
6.8; b. Power-law photon index, fixed at the uniform value for LMXBs found by Irwin et al. (2003). c. In 
units of keV; d. Intrinsic 0.3-10 keV fluxes in units of 10~^^ ergs cm~^ s~^; e. In units of solar, flxed. 
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